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We report on a numerical investigation of the SU(2) gauge theory with gluinos. The low-lying spectrum in 
bosonic and fermionic channels is determined. Improvements of the multi-bosonic algorithm are discussed. 
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1. INTRODUCTION 

Simulating supersymmetric field theory on the 
lattice is not straightforward. As the Poincare 
symmetry is broken on the lattice, so is SUSY. A 
possible way out is to follow the ansatz by Curci 
and Veneziano They proposed to live with 
broken SUSY on the lattice and to reveal the rel- 
evant information in the continuum limit, where 
SUSY should be restored. 

2. ACTION AND ALGORITHMS 

The Curci- Veneziano action oi N = I SYM, 
after integrating out the gluino field, is given by 
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Here the fermion matrix is 



:yv,x 



and Vvu,xii the gauge link variable in the ad- 
joint (triplet) representation. The bare param- 
eters are: the usual gauge coupling (3 — Ag^^ and 
the hopping parameter of gluino K. The factor 
of ^ in front of log det shows that the gluino is a 
Majorana fermion, corresponding effectively to a 
flavour number Nf = i. (For further discussion 
of the lattice action see 1^.) 



'Poster presented by K. Spanderen. 



2.1. Two-step multi-bosonic algorithm 

For the Monte Carlo simulations with dynam- 
ical gluinos we use the two-step variant ^ of the 
multi-bosonic algorithm |^ . This requires smaller 
storage and has shorter autocorrelations. For the 
correction step with the higher order polynomial 
Pn it turned out better to use another method, 
not the one introduced in ref . |^ . In terms of the 
hermitean fermion matrix Q = 75(5 = let us 
define, with n even, the decomposition 
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Using this form, the noise vector 77, necessary in 
the noisy correction step, can be generated from 
the gaussian vector 77' according to 



V = Pn/2{Qy'v' , 

where P,i/2{Q)^^ can be obtained as 
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PnMQr' = 



Pn{Q^) 



P-n{Q^)Pn/2{Q)^ -(5) 



In the last step P^n denotes a polynomial ap- 
proximation for the inverse of P„ on the interval 
[e. A], which covers the spectrum of ~ Q^Q on 
typical gauge configurations. For the calculation 
of the necessary polynomials procedures written 
in Maple are available 
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2.2. Preconditioning and eigenvalue distri- 
butions 

In order to improve the performance of our 
fermion simulation algorithm, preconditioning ac- 
cording to ref. 1^ turned out to be very useful. 
The hermitean fermion matrix is decomposed as 



K = 0.18 j3 = 2.3 



Q = 75*3 



(6) 
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and then we have 

det Q = det Q , with Q = -i^-K^^^MoeM^o -(7) 

The matrix has a smaller condition number 
A/e than . The condition number and its 
fluctuations on different gauge configurations are 
dominated by the minimal eigenvalue. A compar- 
ison of the fluctuations of the lowest eigenvalue of 
and is shown in figure |[ 

Lattice: 4^x8, K = 0.18, /g = 2.3 ^ 
without preconditioning 
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Figure 1. Distribution of the minimal eigenval- 
ues of Q2 and on 4^ • 8 lattice at (/3 = 2.3, 
K = 0.18). 

As a consequence of the smaller condition num- 
bers and smaller orders of the polynomial approx- 
imations the autocorrelations are shorter with 
proconditioning than without it. For an exam- 
ple see tabic ||. 

Random matrix models suggest that the fluc- 
tuations of the minimal eigenvalue are inversely 
proportional to the lattice volume. (For refer- 
ences and a recent summary see ref. |^].) This 
is advantageous for the choice of the interval of 
polynomial approximations [e, A]. Our numerical 
data support the decrease of fluctuations, as is 
shown by figure 0. 
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Figure 2. Distribution of minimal eigenvalues of 
the squared preconditioned hermitean fermion on 
43 • 8 and 6^ • 12 lattice at (/3 = 2.3, K = 0.18). 

3. LOW-LYING SPECTRUM 

A first important question for the numerical 
simulations is to determine the masses of lightest 
states in different channels as a function of the 
gluino mass. These can be compared to the pre- 
dictions of the low energy effective action |P , ||] . 
It is expected that there is confinement as in pure 
gauge theory. In the limit of zero gluino mass 
the states should occur in degenerate supermul- 
tiplets. 

The methods to determine low-lying masses 
can be tested and tuned on "quenched" gauge 
configurations generated with the pure gauge part 
of the action pl,||ll|. 



3.1. Correlations 

The lightest supermultiplet at zero gluino mass 
is presumably a massive chiral multiplet consist- 
ing of a scalar, a pseudoscalar and a spin-^ Majo- 
rana fermion. The bosonic states have analogues 
in QCD and can be made out of gluinos. Since 
the gluinos are in the adjoint representation, let 

Table 1 

Autocorrelation times for the plaquette value on 
43 • 8 lattice at (/? = 2.3, X = 0.18). is the 
deviation norm of the polynomial approximation. 
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Lottice: 6^x-2, )S = 2.3 



Figure 3. The square root of the ratio of the 
connected to disconnected part of the correlation 
for a-eta as a function of 1/K. The hne is a hnear 
extrapolation through the last three points. 

us call them a-eta and a-/o for pseudoscalar and 
scalar, respectively. The states made out of 
gluinos can be generally called gluinoballs. The 
two-fermion correlation functions contain con- 
nected and disconnected parts. The connected 
parts can be considered separately and, in anal- 
ogy with QCD, can be associated with the a- 
pion and a-sigma for pseudoscalar and scalar, re- 
spectively. Near the critical hopping parameter 
Kcr corresponding to zero gluino mass the ratio 
of the connected to disconnected part has to be 
0(1). This is how the dependence of the chi- 
ral symmetry breaking pattern on the number of 
flavours can become manifest. As our numerical 
data obtained with the volume source method [ p^ 
show, for K < Kcr the disconnected part is much 
smaller, but it is increasing for increasing K (see 
figure Extrapolation to a ratio equal to one 
gives a rough estimate Kcr — 0.20. 

3.2. Masses 

Besides the gluinoballs, confined states can also 
be composites of gluons and gluinos. Purely glu- 
onic states are, as usual, called glueballs and the 
spin-i fermionic states made out of gluinos and 
gluons as gluino- glueballs. 

Our first series of Monte Carlo runs with dy- 
namical gluinos has been performed at /3 = 2.3 
and Q.IQ<K < 0.185 on 4^ ■ 8 and 6^ • 12 lattices. 
For the determination of the masses of glueballs 
and gluino-glueballs we used smeared sources [|l3| . 
Preliminary results on the masses of O"*" glueball. 
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Figure 4. Dependence of the lowest bound state 
masses on the hopping parameter X on 6'^ • 12 
lattices at /3 = 2.3. 

a-eta and spin-^ gluino-glueball are shown in fig- 
ure ^. As one can see, in this range of couplings 
the expected degeneracy of bound state masses 
is not yet observed. We are presently extending 
the runs to other parameter values and to larger 
lattices. 
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